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Abstract. In this paper we settle most of the open questions on algorithmic computabil- 
ity of Julia sets. In particular, we present an algorithm for constructing quadratics whose 
Julia sets are uncomputable. We also show that a filled Julia set of a polynomial is always 
computable. 
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1. Foreword 

Computable planar compacts can be visualized on a computer screen with an arbitrarily 
high magnification. Among all computer-generated pictures with mathematical content 
Julia sets of rational mappings occupy, perhaps, the most prominent position. And yet, as 
we have shown in [BY] , some of those sets are uncomputable, and so cannot be visualized. 
In this paper we present an account of the results on computability of Julia sets. We survey 
what was previously known, mostly from our own work, some of it also joint with I. Binder; 
and present new results which settle most of the previously open questions. 

A reader, unfamiliar with the questions of algorithmic computability and complexity, 
particularly as applied to questions in analysis, will find an introduction in §2] An intro- 
duction to the relevant concepts of dynamics of rational mappings awaits in ^ Having 
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thus set the stage, we derive some prehminary results on computabihty of Juha sets in ^ 
In [|5]we present our first new result: 

All filled Julia sets of polynomial mappings are algorithmically computable. 

This answers in the positive a question put to us by John Milnor. In the following [|6]we 
discuss negative results. In particular, we show: 

There exist computable complex parameters c, such that the Julia set of J{z^ + c) is not 
algorithmically computable. 

Previously, we could only show the existence of parameters for which the Julia set is 
uncomputable. Now we can give an algorithm for their construction. The question whether 
such a thing is possible would be invariably asked by our colleagues after our talks on the 
subject. The general feeling was, perhaps, that the uncomputability of Julia sets is in some 
ways connected to lack of computabihty of c's. As we see now, this is not the case. 

In the final ^ we interpret our results, and attempt to describe a toy model for uncom- 
putable Julia sets. We also try to answer the "naive" question: "What would a computer 
really draw, when J{z^ + c) is uncomputable?" Much of the discussion in this section is 
motivated by a problem posed to us by Michael Shub. 

We now have a good understanding of computable properties of rational Julia sets. In 
this paper, however, we completely avoid discussing the related computational complexity 
questions. In this area we know comparatively little. In [ BBY2] . jointly with I. Binder, 
we have shown that computational complexity of quadratic Julia sets can be arbitrarily 
high. Evidently, there is an interplay between dynamical properties of Julia sets, and the 
hardness of drawing their picture on a computer screen. However, as the work [Brv3j of the 
first author suggests, this connection may not be straightforward: it shows that a class of 
maps with "bad" dynamical properties has poly-time computable Julia sets. Undoubtedly, 
many similar surprises still await us in this line of investigation. 

Acknowledgements. It is our pleasure to thank our friend and colleague Ilia Binder 
for the many useful discussions on computabihty of Julia sets. We thank John Milnor 
for posing the question on computabihty of filled Juha sets to us. We are grateful to 
Michael Shub for formulating a question which has motivated much of the discussion on the 
"shape" of uncomputable Julia sets in this paper. We also wish to thank many colleagues 
who have invariably asked us if parameters for uncomputable Julia sets can be produced 
algorithmically. 

2. Introduction to computability 

One of the main goals of computability theory is to classify problems according to 
whether or not they can be solved algorithmically. In fact, such questions existed before 
computers. A famous example is Hilbert's Tenth Problem: 

"Given a diophantine equation with any number of unknown quantities and with rational 
integral numerical coefficients: to devise a process according to which it can be determined 
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by a finite number of operations whether the equation is solvable in rational integers. " 
[Bulletin of the American Mathematical Society 8 (1902), 4 37-479.] 

In other words: 

Is it algorithmically possible to determine if a given diophantine equation is solvable? 

It is fairly clear what an affirmative answer would mean in this method to check 

if an equation has a solution. Giving a negative answer (which turns out to be the correct 
one) requires a more formal definition of "methods" that can be used in the solution - as 
one would need to prove that none of these methods work. A universaly accepted model 
was introduced in 1936 in a seminal work by Turing |Tur] in the form of a Turing Machine. 

2.1. Discrete computability and the Turing Machine. The definition of a Turing 
Machine (TM) is somewhat technical and can be found in all texts on computability eg. 
|Pap[ |Sip| . The computational power of a Turing Machine is equivalent to that of a RAM 
computer, and one can think of it as a program on such a computer. The program can use a 
finite amount of memory at each stage of the computation, but it can always request more, 
and there is no a-priori limit on the amount of memory that the machine uses. In fact, 
there is a general belief, usually referred to as the Church- Turing thesis, which states that 
any computation performed on a physical device can be simulated using a Turing Machine. 
Just as an ordinary computer program, any Turing Machine admits a finite description. 

The definition of a Turing Machines gives a natural way of classifying the computability 
of functions in the discrete setting, such as functions acting on the set of naturals N or the 
set of finite binary strings {0, 1}*. Namely, a function f{x) is computable, if there exists a 
TM which takes input and outputs the value f{x). 

Computable functions are sometimes called recursive. They include simple functions 
such as integer arithmetic operations and lexicographical sorting of strings. They also 
include problems that appear to be difficult in practice, but can be solved nonetheless 
if we are willing to wait sufficiently long. These include, for example, finding the prime 
factorization of an integer and finding the optimal strategy in the game of Go. 

On the other hand, there are many functions that are not computable. One argument to 
see this is a simple counting argument: any TM has a finite description, and hence there 
are countably many TMs. On the other hand, there are uncountably many functions from 
N to N, or even from N to {0, 1} - and thus "most" functions are not computable. It is 
much more iteresting to have specific examples of non-computability. 

One such example is the Halting Problem. The halting function H maps a pair (T, w) 
where T is an encoding of a TM M and w is a, binary input to 1 if the machine M running 
on input w eventually halts, and otherwise. 

Sketch of proof that H is not computable. The proof is by a simple diagonalization argu- 
ment. Suppose there were a TM Mi computing the halting function. Let M2 be the 
following machine: on an input w, M2 uses Mi to compute H{w,w). If H{w,w) = 0, then 
M2 halts, otherwise it goes into an infinite loop. 
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Let W2 be the encoding of M2. What will be the outcome of running M2(w2)? If M2 
halts on W2, then H{w2, W2) = 1, and thus M2 cannot halt on W2 by definition. If M2 fails 
to halt on W2, then H{w2,W2) = 0, and by its definition M2 halts on input W2- In either 
case we arrive at a contradiction. □ 

Consider a predicate A : N x N ^ {0, 1} defined as follows. On an input (x, t), A viewes 
X as an encoding of a pair (M, w) of a TM and an input. A{x, t) is 1 if and only if x gives a 
valid encoding and M halts on w in exactly t steps. It is easy to see that A is a computable 
predicate using a simple simulation. On the other hand, computing the predicate 

B{x) = 3t A{x,t) 

is as difficult as solving the Halting Problem, and thus B is non-computable. This example 
will be useful later on. More generally, a predicate of the form P{x) = 3y R{x, y) for a 
computable predicate R{x, y) is said to be recursively enumerable. Moreover, R{x, y) can 
be modified, so that for every x there exists at most one y such that R{x, y) holds. We 
emphasize this by writing P{x) = 3\y R{x, y) Note that any recursive predicate is also 
recursively enumerable. 

Another explicit example of a non-conputable function is given by the negative solution 
to Hubert's Tenth Problem, which is due to Matiyasevich (see |Matj for details and the 
history of the problem). 

Theorem 2.1. The function that maps an encoding of a diophantine equation E to 1 if E 
is solvable and to otherwise, is non-compuable. 

One of Turing's original motivations for introducing the Turing Machine was classifying 
real numbers into computable and non- computable ones. A number is said to be com- 
putable if there exists a TM that writes its (infinite) decimal expansion digit by digit. An 
equivalent, but slightly less representation dependent is the following definition. 

Definition 2.1. A real number a is said to be computable, if there is a computable function 



: N ^ N such that for all n, 
by Mc- 



a - ^ 



< 2 The set of the computable reals is denoted 



In other words, there exists an algorithm to approximate a with any desired degree of 
precision. As with discrete functions, "most" numbers are non-computable, while most 
"nice" numbers such as vr and e are. It can be shown that Mc with the usual arithmetic 
operations forms a closed real field. 

Here we give an extension of the computable numbers that will be useful later on in the 
paper. 

Definition 2.2. A real number a is said to be right computable, if there is a computable 
function : N — > Q such that 

• the sequence {0(^)} is nonincreasing: 0(1) > 0(2) > . . .; and 

• the sequence {4>{n)} converges to a: lim„^oo0('^) = «• 
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It is obvious that a computable real number is also right- computable. The converse is 
not true in general: 

Proposition 2.2. Right computable numbers form a dense subset inM.\M.c. 

Proof. It is obviously sufficient to present a single right computable number which is not 
computable, as then a dense set can be produced using simple arithmetic manipulations. 
Let P{x) = 3\y R{x,y) be a non-computable predicate on N such that R{x,y) is com- 
putable, as discussed above. Consider the number 



l-X^P(a;)-4- 



a 

x=l 

Then a is non-computable, since computing a would also enable us to compute the pred- 
icate P. On the other hand, a is right computable, as demonstrated by the following 
computable function: 

n n 

</)H = 1 - ^ 5^ y) ■ 4-^ 

x=l y=l 

4>{n) is obviously non-increasing, and 

oo oo oo 

lim (j){n) = 1 - V V R{x, y) ■ 4"^ = 1 - V P{x) ■ 4"^ = a. 

x=l y=l x=l 

□ 

A more detailed discussion on the different extensions of the concept of a computable 
number can be found in [Weij . 

The above definition of computability using Turing Machines directly applies only to 
computability questions for discrete objects. It has to be extended if we want to discuss 
computability of continuous objects such as functions over R or subsets of R^. 

2.2. Oracle computation, computable real functions. The history of defining com- 
putability for real objects probably begins with the work of Banach and Mazur [BMj of 
1937, only one year after Turing's paper. This work has founded the tradition of Com- 
putable Analysis (sometimes also called Constructive Analysis). Interrupted by war, it 
was further developed in the book by Mazur [Maz] . Much research took place in the mid 
1950's in the works of Grzegorczyk [Grzj . Lacombe [Lac], and others. A parallel school 
of Constructive Analysis was founded by A. A. Markov in Russia in the late 1940's. A 
modern treatment of the field can be found in |Kolj and [Weij . 

The definition of computability over the reals presented here falls into this framework. 

Consider the simplest case in which we would like to compute a function / : M — R. On 
an input x, we are trying to compute f{x). As in the case with real numbers, the machine 
M computing / should be able to output f{x) with any given precision 2~". The machine 
M, as well as a practical computer, can only handle a finite amount of information, and 
thus is not capable of reading or storing an entire input x. Instead, it is allowed to request 
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the input x with an arbitrarily high precision. In other words, it has an external tape and 
a command READ(m) which requests a 2 "'"-approximation 0(m) of x to be written on 
this tape. It can then be read by the machine from the external tape. It is convenient to 
take all the approximations from the dyadic set D = : k & Z e N}, as they possess 
a natural finite binary encoding. 

To formally define computability of real functions let us first introduce the notion of an 
oracle: 

Definition 2.3. A dyadic- valued function : N — > D is called an oracle for a real number 
X if it satisfies |0(m) — x\ < 2~'^ for all m. 

An oracle Turing Machine is a TM which can query the value 4>{m) of some oracle (f> for 
an arbitrary m e N. Note that the oracle itself is not a part of the algorithm, but 
rather enters as a parameter. We will use a notation M'^ to empasize the dependence of 
the output of the TM on the values of the oracle. 

To get used to the terminology, imagine a trivial algorithm which given an n e N and a 
good enough approximation of x e IR outputs a 2 ""-approximation of the number 2x. The 
algorithm executes the command 

READ X WITH PRECISION 2-^"+^). 

At this point the user (playing the role of an oracle in the dictionary sense) enters from 
the keyboard a dyadic rational d for which |d — x| < 2"^"+^^. The algorithm proceeds to 
output 2ci as the answer. 

Definition 2.4. Let S he a. subset of R, and let / : 5* — >^ M be a real- valued function on S. 
Then / is said to be computable if there is an oracle Turing Machine M'^{n) such that the 
following holds. If is an oracle for x & S, then for every n G N M'^{n) returns a dyadic 
number q such that \q — f{x)\ < 2~". 

Note that M*^ is supposed to work with any valid oracle 4> for x. The definition generalizes 
trivially to functions with A; > 1 variables. 

Examples of computable functions include most common functions such as an integer 
power, exp(a;), and any trigonometric function. A constant function f{x) = a is computable 
if and only if a is a computable number. 

The oracle terminology allows us to separate the problem of computing the parameter 
X from the problem of computing the function / on a given x. For example, the function 
X I— > a:^ is computable. Hence even if a is a non-computable number, we are still able to 
compute a^, provided we have an oracle access to a. This is despite the fact that is a 
non- computable number. 

A fundamental fact about computable functions in this setting is that computable func- 
tions are continuous: 

Theorem 2.3. Let S <zM}^ , and suppose f : S ^ M. is computable by an oracle machine 
M^. Then f is continuous on S. 
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Proof. Let x G S* and e > be given. Choose an integer m such that 2""^ < 6/2. Let 
(j){n) be an oracle for x such that \(l){n) — x\ < 2^^"+^^ for all n (thus "exceeding" the 
minimum requirement from an oracle). Then M'^{m) otputs a number d G D such that 
\d — f{x)\ < 2"™. It terminates after finitely many steps, and hence is only queried 
up to some finite precision 2~'^. It is now not hard to see that for any x' such that 
\x — x'\ < 2^^~^., there is a valid oracle 0' which agrees with up to precision 2~^. Thus 
for any x' G 5n (x-2-^-^x + 2-*^-i), M'^'(m) outputs the same answer d, and we must 
have \d — f{x')\ < 2^™. Hence for every x' G 5 such that \x — x'\ < 2~^~^, we have 

\f{x) - f{x')\ < \f{x) -d\ + \f{x') -d\< 2— + 2-™ < e. 

□ 

In particular, it shows that discontinuous functions, such as arcsin or xq cannot be com- 
puted by a single machine on the whole domain of definition. 
Same considerations can be used to prove a stronger result: 

Theorem 2.4. In the conditions of Theorem \2. 31 there exists a computable function fi{x, k) : 
X N ^ N such that 

\f{y) - f{x)\ < 2~^ whenever y E S and \y - x\ < 2~^^'^'^\ 

We will refer to the this property by saying that / has a computable local modulus of 
continuity. 

Remark 2.1. In some cases, for example when S = [0, 1], the global modulus of continuity 
(or simply the modulus of continuity) of / on 5 is also computable. That is, we can 
compute a function /i : N — N such that 

(2.1) for any X, y G 5 with \x ~ y\ < 2-"^''^ ^ \f{x)-f{y)\ < 2~K 

More generally, this is true whenever S" is a compact computable set (as will be defined in 
the next section). In particular, this is true whenever S = [a, b] with computable endpoints 
a and b, or when S is the unit circle in M^. 

2.3. Computability of subsets of M'^. Let C M'^ be a compact set. We would like 
to give a definition for K being computable. In the discrete case the distinction between 
computability of functions and sets is not as important, since a set 5* is usually said to be 
computable, or decidable, if and only if its characteristic function xs is computable. The 
same definition would not work over M, since only continuous functions can be computable, 
hence xk would not be computable unless K = (/). 

We say that a TM M computes the set K if it approximates K in the Hausdorff metric. 
Recall that the Hausdorff metric is a metric on compact subsets of M.^ defined by 

dH{X,Y) = inf{e > 0|X C U,{Y) and Y C U,{X)}. 

We approximate K using a class C of sets which is dense in metric dn among compact sets, 
and such that elements of C have a natural binary encoding. Namely C is the set of finite 
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unions of dyadic balls: 



C = <^ y B{di, ri) I where di e ©^ G D 




Members of C can be encoded as binary strings in a natural way. The following definition 
is equivalent to the set comput ability definition given in [Wei] . and in earlier works (e.g. 



Definition 2.5. We say that a compact set if C M'^ is computable, if exists a TM M(m), 
such that on input m, M(m) outputs an encoding of Cm G C such that dni^K, Cm) < 2""^. 

To illustrate the robustness of this definition we present the following two equivalent 
characterizations of computable sets (see e.g. [Brvj ) . The first one relates the definition 
to computer graphics. It is made more precise in the discussion below. The second one 
relates the computability of sets to the computability of functions as per Definition 12.41 

Theorem 2.5. For a compact K dMl' the following are equivalent: 

(1) K is computable as per definition \2.5l 

(2) (in the case k = 2) K can be drawn on a computer screen with arbitrarily high 
resolution, 

(3) the distance function dxix) = inf{|x — y\ | y G K} is computable as per definition 



Let us elaborate further on part (2) of the theorem. A "drawing" P of the set K on the 
computer screen is just a collection of pixels that serve as an accurate description of K (or 
a portion of K, if the image is zoomed- in). We would expect the following properties from 
P: 

• P should include all pixels that intersect with K, this guarantees that we get a 
picture of the entire set P; and 

• P should not include pixels that are "far" from K, for example pixels that are at 
least one pixel diameter away from the set K. 

By switching from the rectangular computer pixels to the mathematically more convenient 
round pixels, we see that to "draw" K one should be able to compute a function fx '■ 
Dx©'=^{0,l} from the family 



Jk then can be used to decide whether to include a round pixel with center d and radius 
r in P. Sample values of the function fx are illustrated on Figure [TJ 



(2.2) 




1 if 5(c/,r)nir^0 

if 5(c/,2 -r) n/sT = 



or 1 otherwise 



2.4. Weakly computable sets. In this section we present a different definition of set- 
comput ability, we call weak computability. If was first introduced by Chou and Ko |CK] . 
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Figure 1. Sample values of the function 

Definition 2.6. We say that a set S is weakly computable if there is an oracle Turing 
Machine M^{n) such that if = 02, • • • , <Pk) represents a point x = {xi, . . . , x^) G K*^, 
then the output of M'I'in) is 

(I lixeK 
(2.3) M'^{n) = l if 5(a;,2-("-i)) n/s: = 

I or 1 otherwise 

Condition (12.31) is similar to condition (12. 2p . The difference is that now we allow x to be 
any point in M*^ (not just D^), and we do not require the machine to output 1 if x is not 
in K but is "close" . It is evident from Figure [2] that Definition 12.61 requires less effort from 
the algorithm computing K than the original definition. Thus, the new definition appears 
to be weaker than the defintion of set computability from last section, but it turns out that 
they are equivalent. 

Theorem 2.6. |Brv2j A compact set K C M.'' is weakly computable if and only if it is 
computable as per definition \2.5i 

It is sometimes easier to use weak computability when proving that a certain set is 
computable. We will use Theorem 12.61 in ^ 

2.5. Set-valued functions and uniformity. The problem of computing Julia sets is 
essentially that of mapping the coefficients of a rational function R{z) to the set Jr. Thus 
we need a notion of computability of set- valued functions to discuss computability questions 
about Julia sets. 

We can now combine the definitions from previous sections to define computability of 
set valued functions. 

Definition 2.7. Let S" be a subset of M.^. Denote by the set of all the compact subsets 
of M^. Let F : S ^ be a set-valued function mapping points in S to compact subsets of 
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Regular 



Weak 




Figure 2. The values of /^^(•,2 ") in the definitions of regular (left) and 
weak set computability 



M^. F is said to be computable on S if there is an oracle TM M'^i' "''^'=(n) that for oracles 
representing a point x = {xi,X2, ■ ■ ■ ,Xk) G S outputs an encoding of a set C„ G C such 
that the Hausdorff distance dH{F{x), C„) < 2~". 

In fact, the computability definitions for real functions, sets, and set-valued functions 
presented above fit in nicely within the much more general framework of Type Two Effi- 
ciency (TTE). See |Wei] . and references therein for more details. In particular. Theorem 
12.31 stating that computable =^ continuous holds in a very broad variety of settings. We 
will only need it in the case of set-valued functions. The proof is very similar to the proof 
of Theorem 12.31 (see e.g. |B Yj ) . 



Theorem 2.7. Suppose F : S C M is computable as per Definition \27' 

continous on S in the Hausdorff metric. 



then F is 



Example. Let the complex plane C be naturally identified with M^. Let d > 1 be an 
integer. Consider the multi-valued function fd=^ : C — C There is no continuous 
single- valued branch of fd on the entire complex plane, hence there is no computable branch 
if fd that is defined on the entire C. There is a computable branch of fd that is defined 
everywhere except for a slit connecting to oo. 

On the other hand, if we view the function fd as a set-valued function that maps a 
number z = r- e^'"'^ to its d roots {r^l'^ ■ e^'^'^l'^, r^l'^ ■ e2'^»(^+i)/«', . . . , r^l'^ ■ e27ri(e+d-i)/d|^ ^j^^^ 
it is not hard to see that fd becomes computable. And indeed, the map : — > is 
continuous in the Hausdorff metric. 
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Note that Definition 12.71 makes sense even when S = {s} is a singleton. In this case 
we say that F is nonuniformly computable on s. Otherwise, we say that F is uniformly 
computable on the set S. 

Remark on the BSS computability model. We note that another approach to com- 
putability of subsets of M'^ has been developed by Blum, Shub, and Smale |BCSSj . It is 
based on the concept of decidability in the Blum-Shub-Smale (BSS) model of real com- 
putation. The BSS model is very different from the Computable Analysis model we use, 
and can be very roughly described as based on computation with infinite-precision real 
arithmetic. Some discussion of the differences between the models may be found in |BrC] 
and |Brv2j . Algebraic in nature, BSS decidability is not well-suited for the study of fractal 
objects, such as Julia sets. It turns out (see Chapter 2.4 of |BCSS ]) that in the BSS model 
all, but the most trivial Julia sets are not decidable. More generally, sets with a fractional 
Hausdorff dimension, including ones with very simple description, such as the Cantor set, 
are BSS-undecidable. 

3. Julia sets of rational mappings 

3.1. Basic properties of Julia sets. An excellent general reference for the material in 
this section is the book of Milnor [Milj . For a rational mapping R of degree deg R = d > 2 
considered as a dynamical system on the Riemann sphere 

R:C^C 

the Julia set is defined as the complement of the set where the dynamics is Lyapunov-stable: 

Definition 3.1. Denote F{R) the set of points z E C having an open neighborhood U{z) 
on which the family of iterates R^\u{z) is equicontinuous. The set F{R) is called the Fatou 
set of R and its complement J{R) = C \ F{R) is the Julia set. 

In the case when the rational mapping is a polynomial 

P(z) = ao + aiz + ■ ■ ■ + adz'^ : C ^ C 

an equivalent way of defining the Julia set is as follows. Obviously, there exists a neigh- 
borhood of cxo on C on which the iterates of P uniformly converge to oo. Denoting A{oo) 
the maximal such domain of attraction of oo we have A{oo) C F{R). We then have 

J(P) = dA{oo). 

The bounded set C \ A{oo) is called the filled Julia set, and denoted K{P); it consists of 
points whose orbits under P remain bounded: 

K{P) = {zet\ sup |P"(^)| < oo}. 

n 

For future reference, let us summarize in a proposition below the main properties of Julia 
sets: 
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Proposition 3.1. Let R : C ^ C be a rational function. Then the following properties 
hold: 

• J{R) is a non-empty compact subset ofC which is completely invariant: R~^{J{R)) = 
J{R); 

• J{R) = J(/?") for all n G N; 

• J{R) has no isolated points; 

• if J{R) has non-empty interior, then it is the whole of C; 

• let U d d be any open set with U fl J{R) 7^ 0. Then there exists n G N such that 
R^{U) D J{R); 

• periodic orbits of R are dense in J{R). 

Let us further comment on the last property. For a periodic point Zq = R^{zq) of period 
p its multiplier is the quantity A = A(zo) = DR^{zq). We may speak of the multiplier of a 
periodic cycle, as it is the same for all points in the cycle by the Chain Rule. In the case 
when |A| 7^ 1, the dynamics in a sufficiently small neighborhood of the cycle is governed by 
the Mean Value Theorem: when |A| < 1, the cycle is attracting {super- attracting if A = 0), 
if |A| > 1 it is repelling. Both in the attracting and repelling cases, the dynamics can be 
locally linearized: 

(3.1) ijiRP{z)) = X-^{z) 

where ip is a conformal mapping of a small neighborhood of zq to a disk around 0. 
In the case when |A| = 1, so that A = e^'"*^, 6* G M, the simplest to study is the parabolic 
case when ^ = n/m G Q, so A is a root of unity. In this case R^ is not locally linearizable; it 
is not hard to see that zq G J{R). The description of the dynamics in a small neighborhood 
of a parabolic orbit will be discussed below in some detail. 

In the complementary situation, two non-vacuous possibilities are considered: Cremer 
case, when R^ is not linearizable, and Siegel case, when it is. In the latter case, the 
linearizing map ip from fl3.ll) conjugates the dynamics of R^ on a neighborhood U{zo) to 
the irrational rotation by angle 6 (the rotation angle) on a disk around the origin. The 
maximal such neighborhood of zq is called a Siegel disk. 

A different kind of a rotation domain may occur only for a non-polynomial rational 
mapping R. A Herman ring A is a conformal image 

u : {z e C\ < r < \z\ < 1} ^ A, 

such that 

RPou{z) = v{e^^'^z), 

for some p G N and G M \ Q. 

The term basin in what follows will describe the set of points whose orbits converge to a 
given periodic orbit under the iteration of R. We will denote Postcrit(-R) the post-critical 
set of R, defined as the closure of the union of the orbits of critical points of R. Fatou 
made the following observation: 
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Proposition 3.2. Let pi, . . . ,pk be a periodic orbit of a rational mapping R. If it is either 
attracting, or parabolic, then its basin contains a critical point of R. 

By a perturbative argument, Fatou then concluded that for a rational mapping R with 
degi? = d > 2 at most finitely many periodic orbits are non-repelling. A sharp bound on 
their number depending on d has been established by Shishikura; it is equal to the number 
of critical points of R counted with multiplicity: 

Fatou-Shishikura Bound. For a rational mapping of degree d the number of the non- 
repelling periodic cycles taken together with the number of cycles of Herman rings is at 
most 2d — 2. For a polynomial of degree d the number of non-repelling periodic cycles in C 
is at most d — 1. 

Therefore, we may refine the last statement of Proposition 13.11 

• repelling periodic orbits are dense in J{R). 
Classical results of Fatou also imply the following: 

Proposition 3.3. Every Cremer point of a rational mapping R as well as every point of 
the boundary of a Siegel disk or a Herman ring is contained in Postcrit(i?). 

By definition, the basin of an attracting or a parabolic point, as well as preimages 
of Siegel disks and Herman rings belong to the Fatou set. Fatou-Sullivan Classification 
Theorem formulated below rules out other possibilities: 

Fatou-Sullivan Classification. For every connected component W C -Fr there exists 
m G N such that the image H = R'^{W) is periodic under the dynamics of R. Moreover, 
each periodic Fatou component H is of one of the following types: 

• a component of the basin of an attracting or a super- attracting periodic orbit; 

• a component of the basin of a parabolic periodic orbit; 

• a Siegel disk; 

• a Herman ring. 

To conclude the discussion of the basic properties of Julia sets, let us consider the sim- 
plest examples of non-linear rational endomorphisms of the Riemann sphere, the quadratic 
polynomials. Every affine conjugacy class of quadratic polynomials has a unique represen- 
tative of the form fc{z) = z"^ + c, the family 

f,{z) = z^ + c, cEC 

is often referred to as the quadratic family. For a quadratic map the structure of the Julia 
set is governed by the behavior of the orbit of the only finite critical point 0. In particular, 
the following dichotomy holds: 

Proposition 3.4. Let K = K{fc) denote the filled Julia set of fc, and J = J{fc) = dK. 
Then: 

• E K implies that K is a connected, compact subset of the plane with connected 
complement; 
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• ^ K implies that K = J is a planar Cantor set. 

The Mandelbrot set C C is defined as the set of parameter values c for which J(/c) is 
connected. 

A rational mapping i? : C ^ C is called hyperbolic if the orbit of every critical point of R 
is either periodic, or converges to an (super-) attracting cycle. The term "hyperbolic" has 
an established meaning in dynamics. Its use in this context is justified by the following 
proposition: 

Proposition 3.5. A rational mapping R of degree d > 2 is hyperbolic if and only if there 
exists a smooth metric /i defined on an open neighborhood of J{R) and constants C > 0, 
A > 1 such that 

> CA" for every z G J{R),n e N. 

As easily follows from Implicit Function Theorem and considerations of local dynamics of 
an attracting orbit, hyperbolicity is an open property in the parameter space of rational 
mappings of degree d>2. 

Considered as a rational mapping of the Riemann sphere, a quadratic polynomial fc{z) 
has two critical points: the origin, and the super- attracting fixed point at oo. In the case 
when c ^ A4, the orbit of the former converges to the latter, and thus fc is hyperbolic. 
Proposition 13.21 implies that whenever has an attracting orbit in C, it is a hyperbolic 
mapping and c E Ai. The following conjecture is central to the field of dynamics in one 
complex variable: 

Conjecture (Density of Hyperbolicity in the Quadratic Family). Hyperbolic pa- 
rameters are dense in M.. 

Fatou-Shishikura Bound implies that a quadratic polynomial has at most one non-repelling 
cycle in the complex plane. Therefore, we will call the polynomial fc (the parameter c, the 
Julia set Jc) Siegel, Cremer, or parabolic when it has an orbit of the corresponding type. 

3.2. Local dynamics of a parabolic orbit. We will describe here briefly the local dy- 
namics of a rational mapping R with a parabolic periodic point p. By replacing R with its 
iterate, if needed, we may assume that R{j)) = p, and R'{p) = 1. The map R then can be 
written as 

R{z) = z + a{z- p)"+^ + 0{{z - p)"+2), for some n G N and a ^ 0. 

Note that the integer n + 1 is the local multiplicity of p as the solution of R{z) = z. 

A complex number i/ G T is called an attracting direction for p if the product au"' < 0, 
and a repelling direction if the same product is positive. For each infinite orbit {R''{z)} 
which converges to the parabolic point, there is one of the n attracting directions u for 
which the unit vectors 

{R''{z) -p)/\R\z) -p\ — > V. 

k—>oo 
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We say in this case that the orbit converges to p in the direction of u. For each attracting 
direction z/, we say that a topological disk U is an attracting petal of R at p if the following 
properties hold: 



• an infinite orbit {R''{z)} is eventually contained in U if and only if it converges to 
p in the direction of u. 

Similarly, U is a repelling petal for R if it is an attracting petal for the local branch of R~^ 
which fixes p. 



Figure 3. A Leau-Fatou flower with three attracting petals (shaded) and 
three repelling petals (emphasized). The attracting and repelling directions 
are also indicated. The arrows show the direction of the orbits in one of the 
petals; the image of this petal is also indicated. 

The petals form a Leau-Fatou Flower at p: 

Theorem 3.6. There exists a collection ofn attracting petals P^, and n repelling petals PJ 
such that the following holds. Any two repelling petals do not intersect, and every repelling 
petal intersects exactly two attracting petals. Similar properties hold for attracting petals. 
The union 



• U 3{p}- 

• i?"(?7) C f/U {p}; 
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The map k conformally transforms the infinite sector between two repeUing directions into 
the plane with the negative real axis removed. In this sector, it changes the map R into 

= t(7 + 1 + 0(1/a/H)) as u; ^ oo. 

Selecting a right half-plane = {Kqz > r} for a sufficiently large r > 0, we have 

ReF{w) > Rew + 1/2, and hence F{H) C H. 

The corresponding attracting petal can then be chosen as the domain k~^{H), using the 
appropriate branch of the inverse. Note, that given the coefficients of the rational mapping 
R, the description of the petal is constructive. Let us formulate this last statement in a 
language suitable for later references: 

Lemma 3.7. For each degree d > 2 there exists an oracle Turing Machine M*^ such that 
the following holds. Let R he a rational mapping of degree d with a parabolic periodic point 
p, with period m and multiplier e^'^**/*. Let n he the number of attracting (and repelling) 
directions at p. The machine Af^ takes as input the values of m, n, s, t and a natural 
number k; it is given oracle access to the coefficients of R and the value of p. It outputs a 
set Lk E C such that the following is true: 

• Lk+i D Lk and U Lk = P is the union of attracting petals of R at p, covering all 
the attracting directions; 

• distj^(Lfc, P) < 2-^ 

The dynamics inside a petal is described by the following: 

Proposition 3.8. Let P be an attracting or repelling petal of R. Then the quotient man- 
ifold P/ z'~^R{z) is conformally isomorphic to the cylinder C/Z. 

In other words, in each of the petals there exists a conformal change of coordinates trans- 
forming R{z) into the unit translation z z + 1. 

Suppose now that the multiplier of the fixed point p is a g-th root of unity, R'{p) = e^'^*^/'', 
where {p, q) = 1. A fixed petal for the iterate R'^ corresponds to a cycle of q petals for R. 
It thus follows that q divides the number n of attracting/repelling directions of p as a fixed 
point of i?^. We make note of the following proposition, due to Fatou: 

Proposition 3.9. Each cycle of attracting petals of a rational mapping R captures an orbit 
of a critical point of R. 

This implies, in particular, that a quadratic polynomial fc with a parabolic periodic point 
C with multiplier 6^"^*^/* has a Leau-Fatou flower at C with a single cycle of q attracting 
petals. 

3.3. Occurence of Siegel disks and Cremer points in the quadratic family. Let us 

discuss in more detail the occurrence of Siegel disks in the quadratic family. For a number 
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G [0,1) denote [ri, r2, . . . , r„, . . .], G N U {oo} its possibly finite continued fraction 
expansion: 

(3.2) [ri,r2,...,r„,...] = !— 



ri 



r2 + 



+ 



r„ + ■ ■ 



Sucli an expansion is defined uniquely if and only if ^ Q. In this case, the rational 
convergents Pn/Qn = [fi,---,fn] are the closest rational approximants of 6 among the 
numbers with denominators not exceeding g„. In fact, setting A = e^'^*^, we have 

- 1| > |A''" - 1| for all < /i< h ^ g^- 

The difference |A^" — 1| lies between 2/g„+i and 27r/g„+i, therefore the rate of growth of 
the denominators g„ describes how well 6 may be approximated with rationals. 

Definition 3.2. The diophantine numbers of order k, denoted 'D{k) is the following class 
of irrationals "badly" approximated by rationals. By definition, 6 G T^lk) if there exists 
c > such that 

Qn+i < cq^~^ 

The numbers can be calculated from the recurrent relation 

g„+i = rn+iqn + Qn-i, with go = 0, gi = 1. 

Therefore, 6 G T>{2) if and only if the sequence {rj} is bounded. Dynamicists call such 
numbers bounded type (number-theorists prefer constant type). An extreme example of a 
number of bounded type is the golden mean 

^, = ^^ = [1,1,1,...]. 

The set 

V{2+) = f]Vk 

k>2 

has full measure in the interval [0, 1). In 1942 Siegel showed: 

Theorem 3.10 ([Sie]). Let R be an analytic map with a periodic point zq E C of period p. 
Suppose the multiplier of the cycle 

A = e^^'^ with e G P(2+), 

then the local linearization equation liS. 1\) holds. 

The strongest known generalization of this result was proved by Brjuno in 1972: 



18 MARK BRAVERMAN, MICHAEL YAMPOLSKY 

Theorem 3.11 ( |Bruj ). Suppose 

(3.3) i?(e) = yi^^^<oo, 

then the conclusion of Siegel's Theorem holds. 



Figure 4. The Julia set of for = [1, 1, 1, 1, . . .] (the inverse golden mean). 

Note that a quadratic polynomial with a fixed Sigel disk with rotation angle 6 after an 
affine change of coordinates can be written as 

(3.4) Pe{z) =z^ + e^^'^z. 

In 1987 Yoccoz |Yocj proved the following converse to Brjuno's Theorem: 

Theorem 3.12 ([Yoc]). Suppose that for 9 G [0, 1) the polynomial Pg has a Siegel point at 
the origin. Then B{9) < oo. 

The numbers satisfying (13.31) are called Brjuno numbers; the set of all Brjuno numbers will 
be denoted B. It is evident that UV{k) C B and thus the set B has full measure in the 
unit circle. On the other hand, it can be shown that its complement is dense-G^. 

The sum of the series (13. 3p is called the Brjuno function. For us a different character- 
ization of B will be more useful. Inductively define 9i = 9 and 9n+i = {l/9n\. In this 
way, 

9n = [^n-, ^n+l-, • • •]• 

We define the Yoccoz' s Brjuno function as 

oo _ 

^{9) = Y,9,92---9n-i\og-. 

n=l 

One can verify that 

B{9) < oo ^ <^{9) < oo. 
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The value of the function $ is related to the size of the Siegel disk in the following way. 

Definition 3.3. Let P{0) be a quadratic polynomial with a Siegel disk Ag 3 0. Consider 
a conformal isomorphism : D i— »• A fixing 0. The conformal radius of the Siegel disk Ag 
is the quantity 

r(0) = |0'(O)|. 
For all other 6 G [0, oo) we set r{6) = 0. 

By the Koebe One-Quarter Theorem of classical complex analysis, the internal radius of 
Ag is at least r{6)/4. Yoccoz |Yoc] has shown that the sum 

$(^) + logr(^) 

is bounded from below independently of 6 E B. Recently, Buff and Cheritat have greatly 
improved this result by showing that: 

Theorem 3.13 ( [BC2] ). The function 9 (^{6) + \ogr{e) extends to M as a 1 -periodic 
continuous function. 

We remark that the following stronger conjecture exists (see |MM Y] ) : 

Marmi-Moussa- Yoccoz Conjecture. [MMY] The function v : 6 ^ $(6') + logr(6') is 
Holder of exponent 1/2. 

Let us remark here, even though we will not use it in the present paper, that in [BY] we 
have demonstrated: 

Theorem 3.14. There exists 9q E B such that the function 6 i— *• $(^) is uncomputahle on 
the domain consisting of a single point {9q} by a Turing Machine with an oracle access to 
9. 

Assuming Marmi-Moussa- Yoccoz Conjecture holds. Theorem 13.141 would be sufficient to 
demonstrate that r{6) is not computable for some values of ^ G T; which in turn, by 
Theorem 13.221 below, would imply non-computability of J{Pg): 

Conditional Implication. // the function 

has a computable modulus of continuity, then it is uniformly computable on the entire 
interval [0, 1]. 

The proof of the above implication uses the following result of Buff and Cheritat ( [BC2] ). 
Lemma 3.15 ( |BC2] ). For any rational point ^ = | G [0, 1] denote, as before, 

Ps{z) = e^-^z + z^, 
and let the Taylor expansion of Pq'^{z) at start with 

P°i(z) = z + Az''+^ + . . . , for q en 
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/ \ 1/9 

Let L{6) = (^) ■ Denote by ^tmnc the modification 0/$ applied to rational numbers 
where the sum is truncated before the infinite term. Then we have the following explicit 
formula for computing v{9): 

log27r 



(3.5) 



v{e) = ^trunc{e) + \ogL{e) + 



Equation (I3.5P allows us to compute the value of v easily at every rational 6 G Qfl [0, 1] with 
an arbitrarily good precision. Assuming that v has a computable modulus of continuity, it 
is computable by a single machine of the interval [0, 1] (see for example Proposition 2.6 in 
|Ko2 j). This implies the Conditional Implication. 
The following conditional result follows: 

Lemma 3.16 (Conditional). Suppose the Conditional Implication holds. Let 9 G [0,1] 
be such that ^{6) is finite. Then there is an oracle Turing Machine Mf computing ^{6) 
with an oracle access to 9 if and only if there is an oracle Turing Machine computing 
r(9) with an oracle access to 9. 

Proof. Suppose that Mf computes ^{9) for some 9. Let M'^ be the machine uniformly 
computing the function v. Then we can use Mf and M'^ to compute log r{9) = v{9) — (^{9) 
with an arbitrarily good precision. We can then use this construction to give a machine 
M2 which computes r{9). 

The opposite direction is proved analogously. □ 




Figure 5. The figure on the left is an attempt to visualize the (uncom- 
putable!) function $, by plotting the heights of exp(— $(6')) over a grid 
of Brjuno irrationals. On the right is the graph of the (conjecturally com- 
putable) function v{x). 

Both figures courtesy of Arnaud Cheritat 
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A note on topological properties of Siegel and Cremer quadratic Julia sets. By 

Proposition 13.31 the Julia set of any Siegel or Cremer quadratic polynomial is connected. 
The following result is due to Sullivan and Douady (see [Sul] ) : 

Theorem 3.17. If the Julia set of a polynomial mapping f is locally connected, then f has 
no Cremer points. Moreover, every cycle of Siegel disks of f contains at least one critical 
point in its boundary. 

Thus, in particular, Cremer quadratic Julia sets are never locally connected. There is a 
vast amount of recent work on pathological properties of Cremer quadratics, and we will 
not attempt to give a survey of results here. Let us only mention a paper of S0rensen |S0r 
in which there is a discussion of the mechanism of non local-connectedness in some such 
sets, which also gives some indication of the visual complexity of pictures of Cremer Julia 
sets. We cannot offer an illustration with a Cremer Julia set to the reader - even though 
we will see that all such sets are computable, no informative pictures of them have been 
produced to this day. 

As for Siegel Julia sets, Petersen |Pet] showed that J{Pe) is locally connected for 9 
of bounded type. A different proof of this was later given by the second author [Yamj . 
Petersen and Zakeri [PZ] further extended this result to a set of angles 6 which has a full 
measure in T. 

On the other hand, Herman in 1986 presented first examples of Pg with a Siegel disk 
whose boundary does not contain any critical points. By Theorem 13.171 the Julia set of 
such a map is not locally-connected. In recent papers of Buff-Cheritat [BCl], and Avila- 
Buff-Cheritat |ABCj it is shown that the boundary dAg of the Siegel disk itself can have 
smoothness just a hair breadth short of analytic in such cases. 

For 9 of bounded type, the boundary dAg is a quasi-fractal Jordan curve (a quasi-circle, 
see |Ahl] ) passing through the critical point cg of Pg. To visualize it, we can use the 
following fact: 

Proposition 3.18. Let 9 be of bounded type, and denote Pn/Qn its continued fraction 
convergents. Let B > be an upper bound on supg„+i/g„. There exist constants K > 0, 
T < 1 which depend only on B , such that 

distH(J^„, dAg) < Kt"", where a„ = {P^icg), 2 = 0,..., g^}. 

The proof is given in the Appendix. Proposition 13.181 gives a recipe for drawing boundaries 
of Siegel disks of bounded type, see, for instance. Figure HI 

Dependence of the conformal radius of a Siegel disk on the parameter. In this 
section we will show that the conformal radius of a Siegel disk varies continuously with the 
Julia set. To that end we will need a preliminary definition: 

Definition 3.4. Let (f/„, m„) be a sequence of topological disks f/„ C C with marked points 
Un G Un. The kernel or Caratheodory conyergence {Un,Un) — > iU,u) means the following: 
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• for any compact K G U and for all n sufficiently large, K G Un', 

• for any open connected set W 3 u, if W G Un for infinitely many n, then W G U. 

The topology on the set of pointed domains which corresponds to the above definition of 
convergence is again called kernel or Caratheodory topology. The meaning of this topology 
is as follows. For a pointed domain {U, u) denote 

(t)^u,u) : O ^ f/ 

the unique conformal isomorphism with 0((7,u)(O) = m, and (0((7,u))'(O) > 0. We again 
denote r{U,u) = |(0(c/,n))'(O)| the conformal radius of U with respect to u. 
By the Riemann Mapping Theorem, the correspondence 

L : {U, u) ^ (f)(u,u) 

establishes a bijection between marked topological disks properly contained in C and uni- 
valent maps : D — >■ C with (f)'{0) > 0. The following theorem is due to Caratheodory, a 
proof may be found in [Pom] : 

Theorem 3.19 (Caratheodory Kernel Theorem). The mapping t is a homeomorphism 
with respect to the Caratheodory topology on domains and the compact-open topology on 
maps. 

Proposition 3.20. The conformal radius of a quadratic Siegel disk varies continuously 
with respect to the Hausdorff distance on Julia sets. 

Proof. To fix the ideas, consider the family Pq with 6 ^ B and denote Aq the Siegel 
disk of Pg. It is easy to see that the Hausdorff convergence J{Pe„) — ^ ^{^0) implies the 
Caratheodory convergence of the pointed domains 

(A,„,0)^(A,0). 

The proposition follows from this and the Caratheodory Kernel Theorem. □ 

In fact, we can state the following quantitative version of the above result. For a pointed 
domain {U,u) denote p{U,u) the inner radius p{U,u) = dist{u,dU). 

Lemma 3.21. Let U be a simply- connected bounded subdomain ofC containing the point 
in the interior. Suppose V G U is a simply- connected subdomain ofU, and dV G B{dU, e). 

Then 

r{U, 0) - r{V, 0) < 4:^/r{U,0)V^. 

Moreover, denote F{x) = 4x/(l + x)^. Then 

The first inequality is based on Koebe Theorem, see e.g. |RZ] for a proof. The left-hand 
side is a standard refinement of Schwarz Lemma. 
An immediate corollary is: 
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Corollary 3.22. Suppose the function r{6) is uncomputable on the set {60} ■ Then the 
function 9 ^ J{Pe) is also uncomputable at the same point. 

Proof. Assume that J{Peo) is computable. Using the output of the TM computing this 
Julia set in an obvious way, for each e > we can obtain a domain V ^ C such that 

V C Ago and dH{dV,dAg^^) < e. 

It is elementary to verify that for every 6 eT, the set J{Pe) C 5(0,2). This implies, by 
Schwarz Lemma, that the conformal radius r{6o) < 2. Hence, by Lemma [3. 2 ![ 

\r{V,0) -r{eo)\ <6 = 8y^. 

Using any constructive version of the Riemann Mapping Theorem (see e.g. |BB] ). we can 
compute r{V, 0) to precision 6, and hence know r{6o) up to an error of 26. Given that 6 
can be made arbitrarily small, we have shown that r(^^o) is computable. 

□ 

We also state for future reference the following proposition: 

Proposition 3.23. Let {6i} be a sequence of Brjuno numbers such that 6i —>■ 6 and 
lim r{6i) = I > 0. Then 9 is also a Brjuno number and r{6) > I. 

Proof. Denote (pi = 0(Ae ,o)- Note that by Schwarz Lemma, the inverse ipi = {(f)i)~^ linearizes 
Pg. on Ag- . By passing to a subsequence we can assure that (pi ^ (p locally uniformly, and 
(p'{0) > /. By continuity, (p~^ is a linearizing coordinate for Pg, so ^ is a Brjuno number. 
Moreover, 0(©) C Ag, and so by Schwarz Lemma r{6) > I. □ 

4. Preliminary results on computability of Julia sets 

4.1. Computability without oracle access to c. It is a natural question to ask how 
easy or how difficult it is to draw a picture of a quadratic Julia set without an oracle access 
to the value of c. As we see below, in such conditions even very simple Julia sets become 
algorithmically uncomputable. Note first the following elementary statement: 

Proposition 4.1. //c G {—00, —2) then fc is hyperbolic, and Jc is a Cantor set. Moreover, 
Jc C -8(0, Pc), where (3c = \/\]'^--c + 1/2 > 2 zs a fixed point of f^,. 

Proof. Let z E C with \z\ = f3c + 6, for some 6 > 0. By the Triangle Inequality, 
\fciz)\ = \z^ + c\> \z^\ + \c\ = \z\^ + c=i(3c + 6f + c> 
> + c + 2(3 J = f^m + 2(3c5 >(3c + A5 . 
It follows immediately that fc{z) 00, and hence Jc C 5(0, (3^. It remains to note that 

c = (3c{l - Pc) < -Pc, and hence /^(c) > Pc- 

□ 

Theorem 4.2. Let c < —2 be an uncomputable real number. Then the Julia set Jc is 
uncomputable by a Turing Machine without oracle access to c. 
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Proof. The fixed point /3c = ■\/l/4 — c + 1/2 of tlie mapping fc is repelling under our 
assumption on c, and lience lies in the Julia set. By the previous proposition, 

Pc = sup \z\. 

Now assume that there exists a Turing Machine M{n) which computes Jc- Use it to 
determine the largest j > such that j ■ is at most 2~"-far from all points in J^. Then 

0<(j-2-"-/5J<2-('^-i), 

hence, jSc is computable. But 

which contradicts the assumption that c is an uncomputable real. 

□ 

4.2. Lack of uniform computability of Julia sets. Another natural question to con- 
sider is whether it is possible to compute all Julia sets, or in particular all quadratic Julia 
sets, with a single oracle Turing Machine M'f'[n). This is ruled out by Theorem 12.71 as the 
dependence c i— J(/c) is discontinuous in the Hausdorff distance. For an excellent survey 
of this problem see the paper of Douady |Dou2] . 

Theorem 4.3 ( |Dou2j ). Denote JI(c) and K(c) the functions c ^ Jc and c ^ re- 
spectively viewed as functions from C to K2 with the latter space equipped with Hausdorff 
distance. Then the following is true: 

(a) if c is Siegel then J(c) is discontinuous at c, but K{c) is continuous at c; 

(b) if c is parabolic then both J(c) and ]K(c) are discontinuous at c; 

(c) if c is neither Siegel, nor parabolic, then both J(c) and K(c) are continuous at c. 

The discontinuity of J at Siegel parameters is not difficult to prove: 

Proposition 4.4. Let & Ai be a parameter value for which fc has a Siegel disk. Then 
the map J(c) is discontinuous at c*. More specifically, let zq be the center of the Siegel 
disk. For each s > there exists c G B{c, s) such that fc has a parabolic periodic point in 
B{zo,s). 

Proof. Denote A the Siegel disk around (0, p its period, and 6 the rotation angle. By the 
Implicit Function Theorem, there exists a holomorphic mapping ( : [/(c*) — ^ C such that 
C(c*) = zo and ({c) is fixed under {fcY- The mapping 

u:c^ DifcYiac)) 

is holomorphic, hence it is either constant or open. If it is constant, all quadratic polyno- 
mials have a Siegel disk. This is not possible: for instance, /1/4 has a parabolic fixed point, 
and thus no other non-repelling cycles. Therefore, u is open, and in particular, there is 
a sequence of parameters c„ — > c=k such that C(cn) has multiplier e^'^*^"/^". Since C(cn) is 
parabolic, it lies in the Julia set of /c„. Hence 

dist^(J(/cJ, J(/eJ) > dist(c„9A)/2 
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for n large enough. □ 

Thus an arbitrarily small change of the multiplier of the Siegel point may lead to an 
implosion of the Siegel disk - its inner radius collapses to zero. 

As an immediate consequence of Proposition 14.41 and Theorem 12.71 we have: 

Proposition 4.5. For any TM M'^{n) with an oracle for c G C denote Sm the set of all 

values of c for which M"^ computes Jc- Then Sm 7^ C. 

In other words, a single algorithm for computing all quadratic Julia sets does not exist. 

4.3. Discontinuity at a parabolic parameter. The discontinuity in JI(c) which occurs 
at parabolic parameter values has found many interesting dynamical implications. The 
proof is very involved, its outline may be found in |Dou2] . It is based on the Douady- 
Lavaurs theory of parabolic implosion. Let us briefly describe its mechanism for the case 
of a quadratic polynomial fc- 

Denote ( a parabolic periodic point of fc with multiplier e^'^'^vln q^yv^ let m G N be 
its period. Let Pa and Pr be an attracting and a repelling petals of fc- Recall that by 
Proposition 13.91 the cycle of images f^"^{PA U Pr), j = 0, . . . , g — 1 forms a full Leau-Fatou 
flower at (. 

By Proposition 13.81 the quotient 

ca = PA/fr - c/z. 

The quotient Ca, is sometimes called the attracting Fatou cylinder. It parametrizes the 
orbits converging under the dynamics of the iterate /™ to the point (. A repelling Fatou 
cylinder Cr ~ C/Z is defined similarly, as the quotient of a repelling petal. 
Let T be any conformal isomorphism Ca — ^ Cr. After uniformization, 

Ca ^ C/Z, Cr ^ C/Z 

t{z) = z + q mod Z for some g G C. Let g^- : Pa ^ Pr be any lift of r; it necessarily 
commutes with /™^. Consider the semigroup G generated by the dynamics of the pair 
(/cA't)- The orbit Gz of a point z G C is independent of the choice of the lift g-r and only 
depends on r. 
Set 

Jic,T) = {z eC such that CzDJc^^}- 
It can be shown that this set is the boundary of 

K(^c,t) = {z & C such that Gz is bounded}. 

Notice that K(^c,t) $ Kc. some of the orbits which converge to ( under fc are thrown into 
the complement (C \ Kc) fl Pr by gr. Holes which thus open in the set Kc motivate the 
use of the term "implosion" . 

The Douady-Lavaurs theory postulates: 

Theorem 4.6. For every r as above and every s > there exists c G B{c, s) such that 
B{Jc, s) D J{c,t)- 
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Figure 6. Before and after a parabolic implosion. The Julia sets (black) 
and filled Julia sets (light gray) of a parabolic quadratic /1/4 (left), and of 
/1/4+e for a small complex e. 

Thus the Julia set of fc grows "bigger" under the perturbation from c to c. 

5. Positive results 
5.1. Computability of filled Julia sets. In this section we show: 

Theorem 5.1. For any polynomial p{z) there is an oracle Turing Machine M^{n) that 
given an oracle access to the coefficients ofp{z), outputs a 2~"'- approximation of the filled 
Julia set Kp = K{p{z)). 

Moreover, 

Theorem 5.2. In the case when p{z) = z^ + c is quadratic, only two oracle machines 
suffice to compute all non-parabolic filled Julia sets: one for c E M., and one for c ^ M.. 

Theorem 15.11 answers in the affirmative the question posed to us by J. Milnor, after we first 
demonstrated the existence of non-computable quadratic Julia sets in |BYj . 
Let us first formulate a general fact: 

Proposition 5.3. Let Q{z) be a complex polynomial. Then there exists a Turing Machine 
M"^ with an oracle input for the coefficients of Q{z) such that the following holds. Consider 
any dyadic ball B = B{x, r) C C, x G D^, r G ©, and let ai, . . . , am be the roots of Q{z) 
contained in B. For any natural number n, the machine M'^ will take n, r, and x as inputs, 
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and will output a finite sequence of complex numbers f3i,. . . ,(3k with dyadic rational real 
and imaginary parts for which: 

• A e 5(x,r + 2-"); 

• each (3i lies at a distance not more than 2""" from some root of Q{z); 

• for every aj there exists Pi with \aj — l3i\ < 2~". 

For a classical reference, see |Wey| ; a review of modern approaches to iterative root finding 
algorithms may be found in [BCSSj . 

For a given polynomial p{z) we construct a machine computing the corresponding filled 
Julia set Kp. We will use some combinatorial information about p in the construction, so 
the algorithm will, in general, vary with the polynomial. Note that all the information we 
will need can be encoded using a finite number of bits. 

• Information that would allow us to compute the non-repelling orbits of the polyno- 
mial with an arbitrary precision, as well as their type: attracting, parabolic, Siegel, 
or Cremer. By Fatou-Shishikura bound, there are at most degp — 1 of them. 

By Proposition 15. 3[ such information could, for example, consist of the list of 
periods fc, of such orbits; and for each i a finite collection of dyadic balls {D- }^!^]^ 
separating the points of the corresponding orbit from the other solutions of the 
equation p'^^ (z) = z. 

• For each (super) attracting periodic orbit ( = {Ci, ■ ■ ■ , Ck}, a finite union of dyadic 
balls = UB{(i, Ti) with the property 

• For each parabolic periodic point with period m and multiplier p/g, the values of 
m, p, q. 

• In the case of a Siegel disc Z), information that would allow us to identify a re- 
pelling periodic point Cd in the same connected component of Kp as D. Again, by 
Proposition 15.31 it is sufficient to know its period, and a small enough dyadic ball 
around it, which separates it from all other points, periodic with the same period. 

5.2. Computing Kp. We are given a dyadic point d G D and an n G N. Our goal is to 
always terminate and output 1 if B{d, 2-'')r]Kp ^ and to output if B{d, 2-2~")nirp = 0. 
We do it by constructing five machines. They are guaranteed to terminate each on a 
different condition, always with a valid answer. Together they cover all the possible cases. 

Lemma 5.4. There are five oracle machines Mext, Mjui, Mattr, Mpar, Mgieg such that 

(1) if d is at distance > | ■ 2^" from Kp, Mext(yd,n) will halt and output 0. If d is at 
distance < 2^" from Kp, Mext{d,n) will never halt; 

(2) if d is at distance < | ■ 2~" from Jp, Mjui{d,n) will halt and output 1. If d is at 
distance > 2 ■ 2~" from Jp, Mjui{d,n) will never halt; 

(3) Mattr(yd,n) halts and outputs 1 if and only if d is inside the basin of an attracting 
orbit of p; 
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(4) Mpar{d,n) halts and outputs 1 if and only if d is inside the basin of a parabolic orbit 
ofp; 

(5) Msieg{d,n) halts and outputs 1 if the orbit of d reaches a Siegel disc, and d is at 
distance > | ■ from Jp. It never halts if d is at distance > 2 ■ 2"" from Kp. 



Proof of Theorem \5.1l given Lemma \5.4\ By Fatou-Sullivan classification it is clear that for 



each {d, n) at least one of the machines halts. Moreover, by the definition of the machines, 
they always output a valid answer whenever they halt. Hence running the machines in 
parallel and returning the output of the first machine to halt gives the algorithm for 
computing Kp. □ 

We now prove Lemma [5. 4[ 

Proof, (of Lemma 15. 4h We give a simple construction for each of the five machines. 

(1) Mext'- Take a large ball B such that p~^{B) d B. Intuitively, we pull the ball back 
under p to get a good approximation of Kp. Let B^ be a 2~*-""'"^^-approximation of 
the set p-^{B). Output iff n B{d, | ■ 2"") = 0. It is not hard to see that this 
algorithm satisfies the conditions on Mext- 

(2) Mjui- By Proposition 15.31 for each k we can compute all periodic orbits of p{z) in 
B with periods j < k, as roots of the equation 

p^{z) — z = 

with an arbitrarily high precision. Moreover, by our assumptions, we have the 
means to distinguish the non-repelling orbits from the repelling ones. 

Let Ck be a finite collection of complex numbers with dyadic rational real and 
imaginary parts which approximate the repelling periodic orbits with periods up to 
k with precision 2^*^"+^\ Output 1 iff d{d,Ck) < ^ ' 2 "- The repelling periodic 
orbits are all in Jp and are dense in this set. Hence the algorithm satisfies the 
conditions on Mjui- 

(3) Mattr- For each attracting orbit ( of period find such that 

i?(p™c(D^-),2-'c) cD^-. 

Let 

1 = 1 + sup and m = Y\ "^f • 

Let Zk be a 2~'-'"'"^)-approximation of p^^^d). If d is inside the basin of an attracting 
orbit C, then Zk will be inside Dq for some k. Output 1 if Zk is inside and at 
least 2~'-far from the boundary of D^. 

(4) Mpar-: We make use of Lemma [3771 Since we can produce arbitrarily good approxi- 
mations of every parabolic periodic point C, of p{z)., we do not need an oracle for the 
value of this point. Let be the sets from Lemma [3.71 corresponding to the point 
C,. Let Zk = p^{d) computed with precision 2~(^+^). We output 1 if Zk is inside 
for some C, and at least 2~^-away from its boundary. 
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(5) Msieg- This is the most interesting case. It is not hard to see that for each k, we 
can compute a union E^. of dyadic balls such that 

Let C* be the center of the Siegel disc (one of the centers, in case of an orbit), and 
let y be the given periodic point in the connected component of C*- We terminate 
and output 1 if Ek separates from y in C (or covers either one of them) for some 
k. 

Clearly, if d is inside the Siegel disc, then the forward images of B{d, | ■ 2~") will 
cover an annulus in the disc that will separate C* from the boundary of the disc, 
and in particular from y. Hence Mgi^g will terminate and output 1. 

On the other hand, if the distance from d to Kp is > 2 ■ 2"", then Ek fl Kp = 
for all k. In particular, Ek cannot separate from y, since they are connected in 
Kp. 

□ 

The proof is simplified in the case of a quadratic polynomial. 

Proof of Theorem \5.2 . . If we assume that p{z) = fc{z) then by the Fatou-Shishikura 
bound, there is at most one non-repelling orbit. By our assumption, it is not parabolic. 
Moreover, if it is a Siegel orbit, then the Julia set is connected. Therefore, any repelling 
periodic orbit will be in the same connected component of Kp as the Siegel disk. 

If c ^ A^, we run Mf.xt and Mj^i- One and only one of them is guaranteed to halt and 
output a correct answer. 

For c G we will use a modified Turing Machine Mgieg- It will compute the set E^ as 
before. If E^ separates the plane, it will use Proposition 15.31 to search for a periodic point 
of period at most k both in the exterior and the interior components oi Ek. If a A; is found 
for which such two orbits are located, or if E^ covers a periodic orbit, it will terminate and 
output 1. 

If c e A^, then we run M^xt, Mj^i, Mhyp, and Msieg. As before, it is easy to see that one 
of them will terminate, and its output will be a correct one. □ 

Corollary 5.5. Denote by V the set of c's for which J^. is parabolic. The function K : 
Kz2+c is continuous in the Hausdorff metric on the set M.\V. 

5.3. Computability of Julia sets in the absense of rotation domains. Similar ideas 
were used in [BBYlj to prove the following theorem: 

Theorem 5.6. Let f be a rational map / : C — C without rotation domains. Then its 
Julia set is computable in the spherical metric by an oracle Turing machine with the 
oracle representing the coefficients of f . The algorithm uses the following non-uniform 
information about each parabolic periodic point C of f with period m and multiplier e^'^'^vln- 
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• a dyadic ball B{w,r) 3 p such that B{w,2r) does not contain any other points 
periodic with period m; 

• the values of m, p, and q. 

Proof. For every natural n we can compute a sequence of rationals {gj} such that 

oo 

(5.1) B{Jf, 2-("+2)) ^ y 2-("+i)) m BiJf, 2-"). 

j=i 

To do that, for each k > n + 2 we compute 2~^-approximations of the periodic points of / 
in C with periods at most k using Proposition [^31 Let M > be some bound on \D'^f"^{z)\ 
in the area of an approximate periodic orbit with period m. Then \Df"^ 
imphes that \Df^{w)\ > 1 for the periodic point w which approximates. In this case we 
add the point to our sequence of rationals. Clearly, for each repelling periodic point of / 
we will eventually obtain in this way a rational point which approximates it with precision 
at least 2"'^""'"^^ Since such points are contained in Jj, and dense there, our sequence has 
the desired property. 

Of course, we can similarly eventually find every attracting orbit ( of f with an arbitrary 
precision. In this case, we will compute a set for this orbit with the same properties as 
before. Set D = U^D^. 

Finally, for each parabolic periodic point C of f let be the sets from Lemma 13.71 Set 
Lk = ^cL^. 

We are now ready to present an algorithm to find a set Cm G C with dist/f(Cm, Jj) < 2~"^. 
Fix m G N. Our algorithm to find Cm € C works as follows. At the k-th step: 

• compute the finite union B^ = ^i=iB{qi,2~^'^~^^^) G C; 

• compute with precision 2~(™"'"^) the complement of the preimage 

f-'^iDULk), 

that is, find Wk G C such that 

dH(W,,C\{f-''{DULk))) < 2-(^'"+3); 

• if Wk C Bk output Cm = Bk and terminate. Otherwise, go to step k + 1. 

By Fatou-Sullivan classification, the algorithm will eventually terminate. Now suppose 
that the algorithm terminates on step k. Since Wk C Bk and J/ C B{Wk, 2"^™+^-') we have 
Jf C B{Cm, 2-(™+3)). On the other hand, U{qi} C J/, and thus Bk = Cm C B{Jf, 2-('"+i)). 

□ 

6. COMPUTABILITY OF JULIA SETS OF SlEGEL QUADRATICS AND NEGATIVE RESULTS 

6.1. Computabilty of r{6) is equivalent to computability of Jg. . Recall the dis- 
cussion of the family Peiz) = z"^ + e'^'^'^^z from §3.31 As before, if ^ is a Brjuno number, 
we denote Aq its Siegel disk, and r{9) the conformal radius of Aq. When 9 ^ B, we set 
r{9) = 0. Let us denote r^up = sup r(6'). We let Jq stand for J{Pe)- 
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Of course, the change in parametrization from c to 6 makes it natural to talk about 
computabilty of Jg by a TM with an oracle for 6, rather than for c. However, these notions 
are obviously equivalent, as c = c{6) is found by the formula: 

(6.1) c = c{e) = A/2 - AV4, where A = e^"*^ 

To address the question of computability of Jg for 9 E B we first make note of the 
following result, proven in |BBYlj : 

Proposition 6.1. Suppose rg is computable by a Turing Machine M'^ with an oracle for 
9. Then so is Jg. 

To clarify the logic of the argument, let us break the proof into two steps: 

Lemma 6.2. Suppose rg is computable by a Turing Machine M'^ with an oracle for 9. 
Then so is the inner raidus p{Ag, 0) = pg. 

Proof. The algorithm works as follows: 

(I) For k E N compute a set G C which is a 2~™-approximation of the preimage 

Pq^{D), for some sufficiently large disk D] 
(II) evaluate the conformal radius r(Dfc,0) with precision 2"^™"'"^) (this can be done, 
for example, by using one of the numerous existing methods for computing the 
Riemann Mapping of a computable domain, see [BB ]): 
(III) as before, denote 

F{x) = 4x/(l + xY, for X G [0, 1]. 

Note that this function is monotone, and let ip{w) = F~^{w). This function is 
computable, and ipi^) = 1- 
Evaluate 

p = ip{rg/r{Dk,0)) 
with precision 2-^""+^^ / p{D,0). If 

|1 -p| < 2-("^+3VpP,0), 

then compute the inner radius p{Dk,0) = around with precision 2"^'""'"^) and 
output this number. Else, increment k and return to step (I). 

Termination. Let K = K{Pg) be the filled Julia set of Pg. Then 

and Dq D Di D D2 D . . .. 

Hence for every 6 > there will be a step k = k{e) after which 

dist{dDk, Jg) <6. 

Since dAg C Jg, by Lemma [3.211 this implies that 

\r{Dk, 0) - r(A,, 0)| = |r(Dfc, 0) - rg\ < A,/r{D,0)V6 ^ 0. 
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Since for every large enough k, the value of 

^(r,/r(Dfc,0)) > 1 -2-(™+4)/p(Z},0), 
the algorithm will eventually terminate on step (III). 

Correctness. Now suppose the algorithm has terminated on step (HI). As C Dk, 
Lemma 13.211 implies that 

1 > — — — > 1 

-p(A,o)- p(AO)' 

and so 

\p{Dk,0)-pe\<2-^^+'\ 

□ 

Lemma 6.3. Suppose pg is computable by a Turing Machine M'^ with an oracle for 9. 
Then so is Jg. 

Proof. The algorithm to produce the 2~"' approximation of the Julia set is the following. 
First, compute a large disk D around with Pe{D) W D. Then, 

(I) compute a set Dk G C which is a 2~*^"+^''-approximation of the preimage Pq^{D)] 
(II) set Wk to be the round disk with radius pe — 2"*^ about the origin. Compute a set 

Bk & C which is a 2~("+^)-approximation of Pg'''{Wk); 
(III) if Dk is contained in a 2^*^"+^)-neighborhood of Bk, then output a 2-("+i)-nei ghborhood 
of Dk \ Bk, and stop. If not, go to step (I). 

A proof of the validity of the algorithm is obvious, and we leave it to the reader. 

□ 




Figure 7. A figure produced by the algorithm of Lemma [631 for 9 = (\/5 — 
l)/2. It does not have the same artistic quality as Figure HI but has a 
guaranteed accuracy up to the selected size of a pixel. 

By Proposition 16. II and Corollary 13.221 we have: 

Theorem 6.4. The conformal radius r{9) is computable by a Turing Machine with an 
oracle for 9 if and only the same is true for the Julia set Jg. 
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Let us make a note: 

Proposition 6.5. Let 9 he of hounded type. Then Jg is computahle hy a TM with an oracle 
for 9. 

Proof. By Proposition I3.18[ pg is computable. The claim follows by Lemma 16. 3[ □ 

Showing that there exist noncomputahle Siegel Julia sets is much more delicate. In the 
next section we will prove that this can happen even if 9 itself is computable. 

6.2. Conformal radius of a Siegel quadratic with a computable 9. The theorem we 
formulate below characterizes the values of r(6') which correspond to computable paiameteis 
9: 

Theorem 6.6. Let r G (0, r^up) he a real numher. Then r = r{9) is the conformal radius 
of a Siegel disc of the Julia set Jg for some computahle numher 9 if and only if r is right- 
computahle. 

Before proving this theorem, let us formulate a corollary: 

Corollary 6.7. There exist computahle values of parameter c, such that the Julia set Jc is 
not computahle hy a TM with an oracle access to c. 

Proof. By Proposition 12.21 there exist right computable numbers r* G [0, rgup] which are 
not in Mc. By Theorem 16.61 r^, = r{9^) for 9^, G Mc- Since 9^, itself is computable, is 
uncomputable by a TM with an oracle access to 9^,. By Theorem 16.41 the Julia set Jg^ is 
uncomputable by a TM with an oracle access to 9^. The claim follows by (16.11) . □ 

Proof of the "only if" direction of Theorem \6.b\ We assume that 9 is computable, and show 
that r{9) is right-computable. Recall that periodic orbits are dense in the Julia set Jg. Let 
Hn be the union of all repelling periodic orbits with periods between 1 and n. By Propo- 
sition 15.31 we can algorithmically find an arbitrarily good approximation of Hn by dyadic 
rationals. 

Jg is connected, and Uif„ is dense in Jg. Thus for every /, there exists rii such that the 
set B{Hj,2~^''~^^^) is connected. Moreover, such rii can be found algorithmically. 

Since Jg separates a from oo, the same is true for B^Hm, 2~^'+^^) provided / is sufficiently 
large. Hence we can compute a strictly increasing sequence {ni}'^i^ C N, and a set Ui <Z C 
with the property 

such that C\Ui has a simply-connected component Wi containing a. 

Using any constructive algorithm for computing the conformal radius |BBj we can ap- 
proximate the k-th term of the sequence 

R, = r{B{W,,2"^'-'^),a) + 5-2'-'^ 
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Hn Je in Hausdorff metric and Uk oo, thus by Lemma [3.211 Rk — * r{9). Moreover, 
{Rk} is a non-increasing sequence. Let pk be a dyadic approximation of Rk that we compute 
so that \pk — Rk\ < 2~^. Let 

rfc = + 3 • 2-\ 
Then {r^} is a computable sequence of dyadic numbers. We have 

hm Tfc = hm pk = hm Rk = r{6), 

fe— >oo fe^oo fe— >oo 

and for each fc, 

rk = Pk + 3-2-'' >Rk + 2-2-'' > Rk+i + 4 • 2-('=+^) > p^+i + 3 ■ 2-('=+^) = r^+i. 

This shows that r{6) is right-computable. 

Note that we know that iJ„ — > is Hausdorff metric, which allows us to conclude 
that Tfc — >■ r(^). However, we do not have (and cannot have) an estimate on the rate of 
convergence of Hn to Jg, and thus cannot obtain an estimate on the rate of convergence of 
fk r{6) and compute r{6). □ 

Proof of the "if" direction of Theorem 16.61 Given a computable sequence {r„} such 
that r„ \ r we claim that we can construct a 6 such that r = r{6). We will be using 
the following three lemmas. The first one is Lemma 3.1 of |BYj . and the second one is 
Lemma 4.2 of |BBY2] . The proofs are outlined in Appendix [Bl 

Lemma 6.8. For any initial segment I = [oq, ai, . . . , a„], write cu = [ao, Oi, . . . , a„, 1, 1, 1, . . . ] . 

Then for any e > 0, there is an m > and an integer N such that if we write (3 = 
[ao, ai, . . . , a„, 1, 1, . . . , 1, N, 1, 1, . . .], where the N is located in the n + m-th position, then 

$(cj) + e < $(/?) < <l>(cu) + 2e. 

Lemma 6.9. Foruo as above, for any e > there is an mo > 0, which can he computed from 
(ao, ai, . . . , a„) and e, such that for any m > niQ, and for any tail I = [a„+m, On+m+i, • • •] 
if we denote 

(3' = [ai , a2, . . . ) a}2, 1, 1, . . . , 1, a^_|_nj, an+m+l; • • -J; 

then 

$(/50 > $(cj) 

Lemma 6.10. Let uo = [ai,a2,...] be a Brjuno number, that is ^{uj) < oo. Denote 
^^k = [O'l, a2, . . . , afc, 1,1,...]. Then for every e > there is an m such that for all k > m, 

^{ujk) < ^{uj) + e. 

Using Proposition 16. 5[ we can get a computable version of Lemmas 16.81 and 16.91 

Lemma 6.11. For any given initial segment I = [ao,ai, . . . , a„] and mo > 0, write u = 
[ao, ai, . . . , a„, 1, 1, 1, . . . ] . Then for any e > 0, we can uniformly compute m > mo, an 
integer t and an integer N such that if we write (3 = [ao, ai, . . . , a„, 1, 1, . . . , 1, N, 1,1,...], 

where the N is located in the n + m-th position, we have 

(6.2) r{uj)-2e <r{(3) <r{uj)-e, 
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(6.3) «l>(/3) > $H, 
and for any 

(6.4) $(7) > ^{u) - 2"". 

Proof. We first show that such m and exist., and then give an algorithm to compute 
them. By Lemma 16.81 we can increase $(to') by any controlled amount by modifying one 
term arbitrarily far in the expansion. 

By Theorem I3.13[ / : 6' 1— > ^{0) + logr(^) extends to a continuous function. Hence for 
any Sq there is a 5 such that |/(x) — f{y)\ < Sq whenever |x — ?/| < 5. In particular, there 
is an mi such that |/(/9) — f{oj)\ < Eq whenever m>mi. 

This means that if we choose m large enough, a controlled increase of $ closely corre- 
sponds to a controlled drop of r by a corresponding amount, hence there are m > mo and 
N such that (16.21) holds. (16.31) is satisfied almost automatically. The only problem is to 
computably find such m and A^. 

To this end, we apply Proposition 16. 5[ Together with Theorem 16. 4[ it implies that for 
any specific m and we can compute r(/3). This means that we can find the suitable m 
and A^, by enumerating all the pairs (m, A^) and exhaustively checking (16.21) and (16. 3p for 
all of them. We know that eventually we will find a pair for which (16. 2p and (16.31) hold. 

Finally, t exists and can be computed by Lemma [6.91 □ 

Lemma 16.101 yields the following lemma. 

Lemma 6.12. The supremum of r(9) over all angles is equal to the supremum over the 
angles whose continued fraction expansion has only finitely many terms that are not 1: 

Tsup = sup r{9). 

e=[ai,a2,...,ak,l,l,...] 

Proof. Let e > be an arbitrary small positive number. By the definition of Vgup there is 
a. 9 = [ai, a2, ■ ■ ■] such that logr(6') > log r^up — e. Denote 

= [(^1,0,2, ■ ■ ■ , Clk, I5 I5 • • •]• 

Lemma [6.1UI states that there is an m such that for k > m, ^{6^) < $(0) + e. Moreover, 
there is a 5 such that whenever 10 — 6*1 < 5 we have 11^(0) — v{6) \ < e. 

6k — >■ 0, hence there is an n > m such that \6n — 0\ < 6. On has the required form, and 
we have 

logr(^„) = viOn) - $(^„) > v{e) - ^{9) -2e = logr(^) -2e> logr,„p - 2,e. 
This shows that we can make r(^„) as close to r^^p as we like. □ 

We are given r = lim \ r„ < Tsup-, hence there is an s and an £ > such that 
'Ts < fsup — 26. By Lemma [6.12[ there is a 70 = [ai, 02, ... , a„, 1, 1, . . .] such that + e/2 < 
^^(70) <rs + e. 
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We are now ready to give an algorithm for computing a rotation number 6 for which 
r{9) = hm \ r„. The algorithm works as follows. On stage k it produces a finite initial 
segment Ik = [oq, . . . , 0^^] such that the following properties are maintained: 

(1) Jo = [ai,a2, . . . ,a„]; 

(2) Jfc has at least k terms, i.e. > k; 

(3) for each k, Ik+i is an extension of 7^; 

(4) for each k, denote 7^. = [4, 1, 1, . . .], then r^+k + 2"(^+^)5 < r{pfk) < Ts+k + '^'^S] 

(5) for each fc, $(7^) > $(7fc-i); 

(6) for each A;, for any extension 

(3 = [h, bm^+i, &mfc+2, • • •], 

> $(7fc) - 2-K 

The first three properties are very easy to assure. The last three are maintained using 
Lemma [6. Ill By this Lemma we can decrease r{'yk-i) by any given amount (possibly in 
more than one step) by extending Ik-i to Ik- Here we use the facts that the r^'s are 
computable and non-increasing. 

Denote 6 = limk^oolk- The continued fraction expansion of 6 is the limit of the initial 
segments Ik- This algorithm gives us at least one term of the continued fraction expansion of 
6 per iteration, hence we would need at most 0{n) iterations to compute 6 with precision 
2~" (in fact, much fewer iterations would suffice). The initial segment 70 can also be 
computed as in the proof of Lemma [6. Ill It remains to prove that, in fact, 6 is the rotation 
number we are looking for. 

Lemma 6.13. The following equalities hold: 

$(0) = lim $(7fc) and r{6) = lim r{'-fk) = r- 

Proof. By the construction, the limit 6 = lim7fc exists. We also know that the sequence 
r{'jk) converges to the number r = lim \ Vk, and that the sequence $(7^) is monotone 
non-decreasing, and hence converges to a value ip {a priori we could have ip = 00). By 
the Caratheodory Kernel Theorem, we have r{6) > r > 0, so ^{6) < 00. On the other 
hand, by the property we have maintained through the construction, we know that ^{0) > 
$(7fc) — for all k- Hence ^{6) > ip- particular, ip < 00. 
From |BC2j we know that 

(6.5) ^ + logr = lim(<l>(7fc) + logr(7fc)) = $(6') + logr(6'). 

Hence we must have ip = ^{0), and r = r{6), which completes the proof. □ 

7. Interpretation of the results 

7.1. How difficult is it to produce a 6 for which Je is uncomputable? As we have 
seen, a value of 6 for which Je is uncomputable can be produced constructively. As we will 
see below, under a reasonable assumption, it is not even hard to do so: 
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Conditional Implication. Assume that the 1-periodic continuous function v : 9 ^ (^{9) + 
\ogr{6) has a computable modulus of continuity h2. this follows, for instance, from 
Marmi-Moussa-Yoccoz Conjecture. Suppose there is a computable sequence ri,r2,... of 
dyadic numbers such that 

• {vi} is non-increasing, ri > r2 > . . ., and 

• limi^ooTi = r. 

Then there is a poly-time computable 9 (and hence a poly-time computable c = c{9)) such 
that r{9) = r. 

Proof. By the assumption, there is a computable function : N ^ N such that 



The proof goes along the lines of the proof of the "if direction of Theorem I6.6[ We 
outline the modifications made to the proof here and leave the details to the reader. The 
key difference is that in the proof of Theorem 16.61 we used Lemma [6. Ill to perform a step in 
decreasing the conformal radius from r(7fc_i) to r{^k)- The algorithm there is basically an 
exhaustive search, which, of course, could take much longer than polynomial time in the 
precision of 7^ to compute. By assuming that v has a computable modulus of continuity, 
we can deal with $(7fc_i) and $(7^) instead of the r(«)'s. We have an explicit formula for 
$ that converges well, and we can compute the continued fractions coefficients to make 
$(7/;) close to whatever we want relatively fast. 

The step of going from 7fc_i to 7^ is as follows. First, we do the following computations: 

• compute dk which is the "drop" in r we are trying to achieve; we want 



• compute using the function a value 5k such that \v{x) — v{ii)\ < dk/8 whenever 
\x-y\< 5k. 

We have no a priori bound on how long these computations would take, but we would still 
like to be computing 9 in polynomial time. To achieve this, we use I's in the continued 
fraction expansion of 9 to "pad" the computation. 

When asked about the value of 9 with precision 2~" which is higher than what the known 
terms of the expansion {Ik-i\ can provide, we do the following: 

• try to compute dk and 5k as above, but run the computation for at most n steps; 

• if the computation does not terminate, output an answer consistent with the initial 
segment [4_i, 1, 1, . . . , Ij; 



• if the computation terminates in less than n steps proceed as described below. 

Note that so far the computation is polynomial in n. For some sufficiently large n the 
computation will terminate in n steps, at which point we will have computed dk and 5k- If 
necessary, we then add more I's to the initial segment to assure that |7fe_i — 7a;| < 5k- 



v{9i) - v{92)\ < 2-" whenever l^i - ^2! < 2"'^^"^ 



4/2 < log(r(7fc_i)) - log(r(7fc)) < 4; 




2n 
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Recall that our goal is to assure that 

4/2 < log(r(7fc_i)) - log(r(7fc)) < 4- 

With the current initial segment for 7^ we have |7fc_i — 7^! < 5^, and hence in the difference 

log(r(7fe_i)) - log(r(7fe)) = $(7^) - $(7^-1) + {v{-fk-i) - v{-fk)) 

the last term is bounded by dk/8. This means that for the current step it suffices to increase 
$(7fc) relative to ^{jk-i) by between |4 and |4- 

Let M be the total length of Ik-i and the I's we have added, and let us extend the 
continued fraction by putting N eN in the M + 1-st term, and all I's further. Increasing 
M if necessary, we can ensure an approximate equality 

$(7fe) ^ $(7fc-i) + a{N) logN 

up to an error of ^4- Let Pm/qm be the M-th convergent of the resulting continued 
fraction. Recall that on an input n we need to compute 9 with precision in time 
polynomial in n. If > l/^/qM, then we do not need to know anything about to 
compute the required approximation. Suppose 2~" < l/^/qM, which means n > loggM/2. 
And we have time polynomial in log qm to perform the computation. 

Note that M — 0{logqM)- It is also not hard to see that a{N) < 2~^/^, so in order 
to have a change by 3dfe/4 we must have A^" > e^(^*^^^), hence by making M sufficiently 
large (depending on the value of dk), we can guarantee that > e^"'^ . This means that 
we can approximate a{N) with the truncated function $ at the M-th convergent of the 
continued fraction. Write Pm/qm = [^i, ^2, • • • , du], and denote 

P — [ai, a2, ■ ■ ■ , Am] • [^2, 03, ... , ajvf] • . . . • [om-i, ^m] • [om]- 

Then /? approximates a{N) within a very small relative error. In particular, we can assure 
that 

/J-(l-^)<«W</3-(l + y. 

In time polynomial in logqM we can compute the exact expression for /3 using rational 
arithmetic: (3 = p/q- Now we can estimate and write it as e^'^'^^^^ in time polynomial in 
log(gM)- From there we can continue by adding enough I's to get Ik and 7^ = [Ik, 1, 1, • • •]• 
By the construction, it would give us the necessary decrease in the value of ^(7^). □ 

7.2. Why is Kq always computable? To provide some intuition why the filled Julia set 
is computable even when the Julia set is not, we propose the following toy model. 
Let A : N ^ {0, 1} be any uncomputable predicate. Consider the set 

( S' U {re^-/^ 1]} for t = (0.444 ■ ■ .)2 e [0, 1) 

ken, dk=i 

y |^g2Wfe| r 1]} for t = 1 

I ifceN, A{k)=i 
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To avoid ambiguity, we always take the finite expansion for dyadic t's. An example of a 
set Qt is depicted on Figure [Hi Firstly, note that if t G (0, 1) is not a computable real, then 
the set Qt is non-computable by a TM without an oracle for t. Moreover, even for a TM 
M'f' equipped with an oracle input for t, the set Qi is clearly non-computable. However, 
when filled, every Qt becomes a computable set - the unit disk. This example suggests 
that filling an uncomputable Julia set we gain computability at the expense of losing the 
long and narrow fjords of the Siegel disk. 




Figure 8. Part of the picture of Qt for t = (0.10101 . . .)2 

7.3. What would a computer really draw? One thing is certain - it is unlikely to 
draw the true picture of an uncomputable Julia set. Indeed, by definition of computability, 
there is no systematic way of distinguishing the picture of Jc from all possible pictures on 
a screen with pixel size 2~". The likelyhood of stumbling upon the correct image when 
the screen size is, say, 800x600 is very remote. The exact nature of the answer would 
depend, of course, on the specific algorithm being used. However, some insight is provided 
by considering the following problem, suggested to us by M. Shub. 
Let 5 be the subset of C x C given by 

Z = {{z,c) : z e Jc}. 

Shub has asked us: 

Is the set Z computable? 

The answer again is "yes" : 

Theorem 7.1. Let d > be any computable real. Then the set 

anc X B{0,d) 

is a computable subset o/C x C. 
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Informally, we may think of projection of^flCx (c — e,c + e) to the first coordinate as 
the picture that a computer could produce when Jc itself is uncomputable. 

To understand how the mechanism of non-computability is destroyed in this case, con- 
sider again the set Qt for t G (0, 1] as the toy model. The set 

W = {{z,t) : 2 e a, t e (0,1]} C C X M 

is computable even though f2j itself is non-computable for t = 1. This happens because in 
the closure of W the "shce" corresponding to t = 1 is 

rh 

feGN 

This set "masks" the computational hardness of Qi, and makes W computable. 

Computability of the set ^. We prove Theorem 17.11 by showing that ^ is weakly com- 
putable (Definition 12.61) . 

We will need the following lemma. 

Lemma 7.2. For any point {z, c) in the complement of the closure z converges to an 
attracting periodic orbit of fc : z ^ z"^ + c. 

The proof of the lemma will require us to recall the nature of discontinuities in the 
function J(c), particularly the theory of parabolic implosion. We will give it in the end of 
the section. 

The following lemma allows us to "cover" all points that belong to ^. 

Lemma 7.3. There is an algorithm Ai{n) that on input n outputs a sequence of dyadic 
points pi,p2, ■ ■ ■ E C X <C such that 

oo 

5(5,2-("+3)) ^ |j5(p^.,2-("+2)) c 5(a2-("+i)). 
i=i 

Proof. The repelling periodic orbits of fc are dense in J^- Hence, the set 

Srep = {{z,c) : z is in a repelling periodic orbit of fc} 

is dense in Srep is a union of a countable number of algebraic curves given by the 
constraints 

r f^{z) = z 
I l(/r)'Wl>i 

For each m we can compute a finite number of points p"^, . . . approximating such 
that 

5(^,^p,2-("+=^)) C |j5(p™,2-("+2)) ^ B(S™p,2-("+i)). 
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We have _ 

o 

-li = Srep — \^ S,^ 



m 

rep' 

m=l 



Hence the computable sequence p\, ■ ■ ■ ^vl^iVii ■ ■ ■ )Pr2' • • • ^pT^ ■ ■ ■ yPrln^ ■ ■ ■ satisfies the con- 
ditions of the lemma. □ 

Corollary 7.4. There is an oracle machine Mf^''^^{n), where 0i is an oracle for z ^ C 
and (j)2 is an oracle for c G C, such that Mf^''^'^ always halts whenever d{{z,c),Z) < 2~^"'~^'^^ 
and never halts if d{{z,c),Z) > 2~". 

Proof. Query the oracles for a point p G C x C such that d{p, {z, c)) < 2~^^~^^\ Then run 
the following loop: 

i ^ 
do 

i^ i+1 

generate Pi using Ai{n) from Lemma [7.31 
while d{p,p,) > 2-("+2) 

If d((z,c),a) < 2-("+^), then d{p,Z) < 2-("+3), hence by Lemma O there is an i such 
that d{p,pj) < 2"'^"'+^), and the loop terminates. If d{{z,c),Z) > 2"", then d{p,Z) > 
2-n _ 2-(n-4) ^ ^ 5 . 2-{n+i). Reuce, by Lemma EH P i 2-("+i)) for all i, and the 
loop will never terminate. □ 

The following lemma allows us to exclude points outside ^ from ^. 

Lemma 7.5. There is an oracle machine M^^'''^^, where 0i is an oracle for z G C and 02 
is an oracle for c G C, such that M^^'*^^ halts if and only if z converges to an attracting 
periodic orbit (or to oo) under fc'-z^z^^c. 

Proof. M2 is systematically looking for an attracting cycle of fc. It also iterates fc on z 
with increasing precision and for increasingly many steps until we are sure that either one 
of the two things holds: 

(1) the orbit of z converges to 00; or 

(2) we find an attracting orbit of fc and the orbit of z converges to it. 

If the search is done systematically, the machine will eventually halt if one of the possi- 
bilities above holds. It obviously won't halt if neither holds. □ 

Proof of Theorem 7.1. The algorithm is: Run the machines Mf^''^'^{n) from Corollary 



17.41 and M2 from Lemma (7.51 in parallel. Output 1 if Mi terminates first and if M2 
terminates first. 

First we observe that Mi(n) only halts on points that are 2~"-close to 5, in which case 
1 is a valid answer according to Definition 12.61 Similarly, M2 only halts on points that are 
outside in which case is a valid answer. Hence if the algorithm terminates, it outputs 
a valid answer. It remains to see that it does always terminate. Consider two cases. 
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Case 1: {z,c) E In this case d{{z,c),Z) = < 2^("+^), and the first machine is 
guaranteed to halt. 

Case 2: {z,c) ^ By Lemma \7.2\ z converges to an attracting periodic orbit of /c in 
this case, and hence the second machine is guaranteed to hah. □ 

Proof of Lemma 17.21 Suppose z ^ Jc and the orbit of z does not belong to an attracting 
basin. By the Fatou-SuUivan classification (see e.g. [Milj ) . there exists G N such that 
w = fc{z) belongs to a Siegel disk or to the immediate basin of a parabolic orbit. Our 
aim is to show that for an arbitrary small 5 > 0, there exists a pair (5, c) G C x C with 
l^; — 5| < S, \c — c\ < 6, and for which z G Jc- We will treat the Siegel case first. 

The case when w lies in a Siegel disk. Let us denote A the Siegel disk containing w, and 
let m G N be its period, that is, the mapping 

: A ^ A 

is conjugated by a conformal change of coordinates : A ^ D to an irrational rotation of 
D. 

By Proposition 14.41 we have the following. Denote ( = 0^^(O) G A the center of the 
Siegel disk. For each s > there exists c G -B(c, s) such that fc has a parabolic periodic 
point ( of period m in B{(, s). In particular, Jg is connected, and B{(, s) H Jc 0. 

Consider now the /^-invariant analytic circle 

Sr = (l)-\{z = re^^'\ e G [0,27r)}) 
which contains w. Let e > be such that 

i?(^,e)c/,^(i?(z,5))nA. 
Set B = B{w,e/2) and let n G N be such that the union 

U fr{B) D Sr. 

0<i<n 

By Proposition 14.41 for all 5 > small enough, there exist c G B(c,S) for which is 
connected and there is a point of Jg inside the domain bounded by Sr. Since repelling 
periodic orbits of fc are dense in dA, again for 6 small enough, there are points of on 
the outside of Sr as well, and so there exists a point ^ E Jc ri Sr. By construction, there 
exists j G N such that fi{B{z, 5)) 3 ^. By invariance of Julia set, if c is close enough to c 
we have B{z, 5) fl Jc 7^ 0, and the proof is complete. 

The case when w lies in a parabolic basin. Denote ( the parabolic periodic point of fc 
whose immediate basin contains w, and let m G N be its period. We employ the notations 
of S 

Recall, that by Theorem 14.61 for every s > there exists c G B{c, s) such that B{Jc, s) D 

J(c,t)- 
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Since ( G Jc, and Jc is connected, there exists a point m G Jc H Pr. Let w G C^i be the 
orbit of w, and let u G C/j be the orbit of u. Choose the translation r : Ca Cr so that 
t{w) = u. Then J(c,t) 3 2;, and the claim follows by Theorem 14.61 

Appendix A. Proof of Proposition 13.181 

Siegel quadratic Juha of bounded type sets may be constructed by means of quasicon- 
formal surgery (cf. [Doul] ) on a Blaschke product 

Uz) = e2--W^2^^_ 

This map homeomorphically maps the unit circle T onto itself with a single (cubic) critical 
point at 1. The angle r(7) can be uniquely selected in such a way that the rotation number 
of the restriction ^(/^It) = 7- 
For each n, the points 

form the n-th dynamical partition of the unit circle. The following result is due to Swiatek 
and Herman (for the proof see e.g. Theorem 3.1 of }dFdMj ): 

Theorem A.l (Universal real a priori bound). There exists an explicit constant B > 1 
independent of'-f and n such that the following holds. Let 7 G M \ Q and n G N. Then any 
two adjacent intervals I and J of the n-th dynamical partition of f^ are B -commensurable: 

< 1^1 < B\I\- 

Proposition A. 2 ([He]). For each hounded type 7 = [oq, . . . , Ofc, . . .] the Blaschke prod- 
uct /-y is Ki-quasisymmetrically conjugate to the rotation : x ^ x + 7modZ. The 

quasisymmetric constant may be taken as Ki = (2 sup Oj)"'^'^'^ . 

Let us now consider the mapping \E' which identifies the critical orbits of and P-^ by 

^■.f;{i)^p;{c,). 

We have the following (see, for example. Theorem 3.10 of [YZj ) : 

Theorem A. 3 (Douady, Ghys, Herman, Shishikura). The mapping \& extends to a 
K-quasiconformal homeomorphism of the plane C which maps the unit disk D onto the 
Siegel disk A^. The constant K may he taken as the quasiconformal dilatation of any 
global quasiconformal extension of the K^-qs conjugacy of Proposition In particular, 
K < 2Ki. 

Elementary combinatorics implies that each interval of the n-th dynamical partition 
contains at least two intervals of the {n + 2)-nd dynamical partition. This in conjunction 
with Theorem I A . 1 1 implies that the size of an interval of the (n + 2)-nd dynamical partition 
of is at most r" where 
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Hence, setting 

= {P^{Cj), i = 0, ...,g„+2}, 

by Theorem \A.3\ 



Appendix B. Proof of Lemmas I6.8L 16.91 and 16.101 

We present outlines of proofs for Lemmas 16.81 16.91 and I6.10[ The complete proofs of the 
intermediate lemmas can be found in |BBY2j and [BYj. For convenience, we restate the 
lemmas here: 

Lemma 16.81 For any initial segment / = [oq, Oi, . . . , a„], write u = [oq, ai, . . . , a„, 
1, 1, 1, ... ]. Then for any e > 0, there is an m > and an integer N such that if we write 
P = [ao, ai, . . . , a„, 1, 1, . . . , 1, A^, 1, 1, . . .], where the is located in the n + m-th position, 
then 

Lemma 16.91 For u as above, for any e > there is an mo > 0, which can be 
computed from ( ) and £, such that for any m > mo, and for any tail / = 

[ttn+m, fln+m+l, • • •] if WC dcUOtC 

(3^ = [ai, a2, . . . , fln, 1, 1, . . . , 1, Cln+m-i Q-n+m+li ■ ■ \i 

then 

Lemma 16.101 Let uj = [ai,a2,...] be a Brjuno number, that is ^{uj) < oo. Denote 
= [^1, 02, . . . , ttfc, 1,1,...]. Then for every e > there is an m such that for all k > m, 

Denote 

1 

$ (iu) = ^{iu) - ao{uj)ai{u) . ..an+m-i{^) log- 



The value of the integer m > is yet to be determined. Denote 

= (ao, ai, . . . , a„, 1, 1, . . . , 1, N, 1,1,...). 

The following estimates are proven by induction. 

Lemma B.l. For any N, the following holds: 
(1) For i < n + m we have 
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(2) for i < n + m, 

(3) for i < n + m, 

(4) for i < n + m — 1, 



log 



log 



log; 



^ 2i—{n+m). 



log 



< 



The estimates yield the following. 



Lemma B.2. For any uj of the form as in lemma lUTSl and for any e > 0, there is an mo > 
such that for any N and any m > mo, 

Proof. (Sketch). The Yl the expression for $(/?^) converges, hence there is an mi > 1 
such that the tail of the sum J2i>n+mi o:i<^2 ■ ■ ■ oii-i log ^ < ^- K can be shown that 

• for a sufficiently large mo > mi, if m > mo, then for any N the influence on the 
sum of the "head" elements is very small: 

n+mi — 1 ^ n+mi — 1 ^ 

5^ ai(n...a,_i(/?^)log^-^- «i(/?^)...a,„i(/5^)log^ 

• for the "tail" terms, for i > n + mi such that i n + m, 

ai(/3^)...a,_i(/3^)log ' 



e 

T6- 



< e. 



. . .a,_i(/?i)log^^ 

After the change each term of the tail could increase by a factor of e at most. The value 
of the "tail" starts at the interval (0, ^], hence it remains in the interval (0, If], and the 
change in the tail is bounded by < f|- 

So the total change in $~ is bounded by 

change in the "he^J" + change in the "tail" < ^ + ^ = ^ 

^ ^ 16 16 4 



□ 



Lemma IB. 21 immediately yields: 



Lemma B.3. For any e and for the same mo{e) as in lemma W7^ for any m > mo and 
N, 
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Denote ^^(u) = ao(uj)ai(uj) . . .an+m-i{^) log — ^-rr = ^{^) — Using the esti- 

mates IB. II one can prove the following: 

Lemma B.4. For sufficiently large m, for any N , 

$i(/3^+^)-<l>i(/3^)<|. 

Since $ = $^ + summing the inequalities in Lemmas IB.3I and IB.4I yields: 
Lemma B.5. For sufficiently large m, for any N , 

- <l'(/3^) < e. 

It is immediate from the formula of that: 
Lemma B.6. 

lim m3^) = oo. 

N^oo 



We are now ready to prove Lemma 16. 8[ 

Proof of Lemma \6.8i . Choose m large enough for lemma IB. 51 to hold. Increase by one 
at a time starting with = 1. We know that = ^{uj) < ^{uj) + e, and by Lemma 

IB.6i there exists an M with ^{P^'^) > (^{cu) + e. Let be the smallest such M. Then 
$(/?^"^) < ^i^) + e, and by lemma [R5] 

<l'(/3^) < <l'(/5^~^) +e< <^{uj) + 2e. 

Hence 

<l>(c<;) +e< ^{{3^) < $(cj) + 2e. 
Choosing (3 = (3^ completes the proof. □ 

The second part of the following Lemma follows by the same argument as Lemma IB. 21 by 
taking iV > 1 to be an arbitrary real number, not necessairily an integer. The first part is 
obvious, since the tail of uj has only I's. 

Lemma B.7. For an uj = (3^ as above, for any e > there is an > 0, such that for 
any m> mo, and for any tail I = [a„+m, o,n+m+i, ■ ■ ■] if we denote 

13' = [ai 



then 



and 



ai(/3>2(/5^)...«.-i(/5^)log^^< 



i>n+m 



n+m— 1 

E 

i=l 



• • • log — ^ - . . . lo, 



< e. 



We can now prove Lemma 16.91 
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Proof of Lemma \6.9[ Applying lemma [R7l with | instead of e, we get 

- $(cj) = ^{"liead"(/?0 - "head"(cj)} + ^{"tail"(/50 - "tail"(cu)} > 
-- - ViWH} >---- = -e. 

Sketch of the proof of Lemma \6.1(A We divide the sum for $(ti;), 

s ^ oo ^ 

'^(^^) = y]"i(^) • • •«i-i(^)log^rT+ ai{uj) . . .ai^i{uj)\og——-, 



□ 



^ V ' ^ ^ 



"head" "tail" 
so that "tail" < e/16. Using the estimates from Lemma [B. II one can show that modifying 
u) to uJk for some appropriately chosen k ^ s will satisfy: 



1 1 

^^^1(0;^) . . log — -. — r < ai{uj) . . . log — -— + e/16, 

^ ai{ujk) ^ ai{uj) 

1=1 ^ ' 1=1 ^ ' 

since the relative error in the "head" terms can be made arbitrarily small; 
• for i > s, 

1 1 
• • • (^i-i{^k) log — r < 9 ■ . . . log —-- + T^l'^ . 

Note that for %> k the 2~^/^ term alone dominates the expression on the left. 
Finally for above, 

^(cufc) < "head"(cufc) + "taiF'K) < "head"(w) + e/16 + 9 "tail" (w) + 2"^+^ < 

"head" (cu) + e/16 + 9^/16 + 2""+^ < $(0;) + e, 
for a sufficiently large s. 



□ 
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